Vortex Formation by Interference of Multiple Trapped Bose-Einstein Condensates 
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We report observations of vortex formation as a result of merging together multiple 87 Rb Bose- 
Einstein condensates (BECs) in a confining potential. In this experiment, a trapping potential is 
partitioned into three sections by a barrier, enabling the simultaneous formation of three indepen- 
dent, uncorrelated condensates. The three condensates then merge together into one BEC, either 
by removal of the barrier, or during the final stages of evaporative cooling if the barrier energy is low 
enough; both processes can naturally produce vortices within the trapped BEC. We interpret the 
vortex formation mechanism as originating in interference between the initially independent conden- 
sates, with indeterminate relative phases between the three initial condensates and the condensate 
merging rate playing critical roles in the probability of observing vortices in the final, single BEC. 

PACS numbers: 03.75.Lm, 03.75.Kk, 03.65.Vf 
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In a superfluid, long-range quantum phase coherence 
regulates the dynamics of quantized vortices Q, Q and 
provides routes to vortex formation that are inaccessible 
with classical fluids. For example, in dilute-gas Bose- 
Einstein condensates (BECs), quantized vortices can be 
created using quantum phase manipulation 0, 0] . Vor- 
tices in BECs have also been created using methods more 
analogous to those in classical fluid dynamics la , namely 
through rotating traps E l [5, IS El, turbulence 0], and 
dynamical instabilities |lllll2| . Yet in contrast with clas- 
sical fluid dynamics, to our knowledge vortex generation 
via the mixing of initially isolated superfluids remains 
an unexplored research area. Due to the availability and 
relative ease of microscopic manipulation and detection 
techniques, BECs are well-suited to address open ques- 
tions regarding superfluid mixing and associated vortex 
generation, along with the possible accompanying roles 
of phase-coherence and matter-wave interference. 

In this Letter, we describe experiments demonstrating 
that the mixing or merging together of multiple conden- 
sates in a trap can indeed lead to the formation of poten- 
tially long-lived quantized vortices in the resulting BEC. 
We ascribe the vortex generation mechanism to matter- 
wave interference between the initially spatially isolated 
but otherwise identical BECs, and show that vortex for- 
mation may be induced even for slow mixing time scales. 
While it is now well-known that matter-wave interference 
may occur between BECs |l3j |. our experiment demon- 
strates a physical link between interference and vortex 
generation, providing a new paradigm for vortex forma- 
tion in superfluids. We emphasize that no stirring or 
phase engineering steps are involved in our work, nor 
are any other means for controllably nucleating vortices 
in the trapped atomic gas; the vortex formation process 
itself is stochastic and uncontrollable, and depends on 
relative condensate phases that are indeterminate prior 
to condensate mixing. The vortex formation mechanism 
identified here may be particularly relevant when defects 
or roughness are present in a trapping potential, or when 



multiple condensates are otherwise joined together. Our 
experiment may also illuminate aspects of vortex forma- 
tion at site defects in other superfluids, for which micro- 
scopic studies may be exceedingly difficult and questions 
regarding vortex formation mechanisms are unresolved. 

To illustrate the basic concept underlying our experi- 
ment, we first consider our atom trap, which is formed by 
the addition of a time-averaged orbiting potential (TOP) 
trap [lij l and a central repulsive barrier created with blue- 
detuned laser light that is shaped to segment the har- 
monic oscillator potential well into three local potential 
minima. Figure ^a) shows an example of potential en- 
ergy contours in a horizontal slice through the center of 
our trap. We will assume throughout the ensuing descrip- 
tions that the energy of the central barrier is low enough 
that it has negligible effect on the thermal atom cloud; 
such is the case in our experiment. However, the central 
barrier does provide enough potential energy for an inde- 
pendent condensate to begin forming in each of the three 
local potential minima from the one thermal cloud. If 
the central barrier is weak enough, condensates with re- 
pulsive interatomic interactions will grow together during 
evaporative cooling; if the barrier is strong enough, the 
condensates will remain independent. In this latter case, 
the central barrier height may be lowered while keeping 
the condensed atoms held in the TOP trap. Overlap and 
interference between the heretofore independent conden- 
sates would then be established as the condensates merge 
together into one. We have examined both scenarios. 

Depending on the relative phases of the three interact- 
ing condensates and the rate at which the condensates 
merge together (via either process), the final merged 
BEC may have nonzero net angular momentum about 
the vertical trap axis, as we now describe. We first recall 
that the relative phase between two independent super- 
fluids is indeterminate until an interference measurement 
is made. However, when interference occurs, a directional 
mass current will be established between the superfluids. 
A relative phase can then be determined, but it will vary 
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FIG. 1: (a) An illustration of the potential energy profile of 
our three-well trap in a horizontal slice through the trap cen- 
ter. This contour plot represents the addition of our harmonic 
trap with the measured intensity profile of our optical barrier 
beam, scaled to a potential energy, (b) A computer-generated 
profile of the binary mask used to create our optical barrier, 
where white represents the transmitting area, (c) An image of 
our optical barrier. (d),(e) Phase-contrast images of trapped 
condensates representing integrated optical depths along the 
vertical trap axis. Each shows an area of 85 /im per side, 
as do (a) and (c). In (d), three independent condensates are 
generated in the presence of a strong barrier beam with 170 
fj,W. (e) With 45 /iW in the barrier beam, the initial conden- 
sates have merged together during evaporative cooling. Note 
that even this weak barrier displaces atoms from the center 
of the BEC, so the observable density dip in this image does 
not signify the presence of a vortex core. 

randomly upon repeated realizations of the experiment 
|l5lll6| . In our experiment, when the initial condensates 
merge together, fluid flow in the intervening overlap re- 
gions is established; a straightforward model of the mass 
current for two overlapping but otherwise uncorrelated 
states shows that the direction of fluid flow depends on 
the sine of the relative phase between the states . 

When our three condensates are gradually merged to- 
gether while remaining trapped, fluid flow that is simul- 
taneously either clockwise or counter-clockwise across all 
three barrier arms may occur with finite probability. For 
ease of this discussion, and keeping in mind that only 
relative phases carry physical meaning, we imagine that 
the condensates formed in the three local minima of Fig. 
^a) can be labeled with phases <f>j, where the indices 
j=l, 2, and 3 identify the condensates in a clockwise 
order, respectively. Upon merging of the three conden- 
sates, if it turns out that for example, 02 — 01 = 0.77T 
and 03 — 02 = 0.87T (thus necessarily 0i — 3 = 0.57T since 
each <j)j must be single valued), then some finite amount 
of clockwise-directed fluid flow will be established be- 
tween each pair, hence also for the entire fluid. More 
generally, if the three merging condensates have relative 
phases 02 — 0i, 03 — 02, and 0i — 03 that are each simul- 
taneously between and 7T, or each between 7r and 2-7T, 
the resulting BEC will have nonzero angular momentum, 
which will be manifest as a vortex within the BEC. By 
examining the full range of phase difference possibilities, 
the total probability P v for a net fluid flow to be estab- 
lished in either azimuthal direction is determined to be 
P v = 0.25, given statistically random phase differences 
for each experimental run. P v is thus the probability for 
a vortex to form as the three condensates merge together. 

Absent from the above description is an analysis of 
the phase gradients of the three condensates during the 



merging process, which will lead to transient interference 
fringes in overlapping condensates. These fringes may 
decay to numerous vortices and antivortices, and possibly 
vortex rin gs, similar to instabilities of dark solitons in 

becs miiiEiii 

For rapidly merged condensates, 
we may thus expect the observation of multiple vortex 
cores in an BEC image, or to find a value of P v greater 
than 0.25. Yet as the condensates are merged together 
more and more slowly, the dependence of P v on phase 
gradients becomes negligible and P v should approach a 
limiting value of 0.25. 

Our experiment is designed to study the presence of 
vortices in an \F—1, mp=—l) 87 Rb BEC subsequent to 
the merging of three condensates such created in a three- 
segment potential well. To create just a single BEC in 
a trap without a central barrier, we first cool a thermal 
gas to just above the BEC critical temperature in an ax- 
ially symmetric TOP trap with radial (horizontal) and 
axial (vertical) trapping frequencies of 40 Hz and 110 
Hz, respectively. We then ramp the TOP trap magnetic 
fields such that the final trap oscillation frequencies are 
7.4 Hz (radially) and 14.1 Hz (axially). A final 10-second 
stage of radio-frequency (RE) evaporative cooling pro- 
duces condensates of ~4xl0 5 atoms, with condensate 
fractions near 65% and thermal cloud temperatures of 
^22 nK. The chemical potential of such a BEC is fcsx8 
nK, where ks is Boltzmann's constant. 

To create instead three isolated condensates in a seg- 
mented trap, we modify the above procedure by ramping 
on an optical barrier immediately before the final 10-s 
stage of evaporative cooling in the weak TOP trap. The 
barrier itself is formed by illuminating a binary mask, il- 
lustrated in Fig.njb), with a focused blue-detuned Gaus- 
sian laser beam of wavelength 660 nm. After passing 
through the mask, the beam enters our vacuum chamber 
along the vertical trap axis. The mask is imaged onto 
the center of the atom cloud with a single lens. Due to 
diffraction, the beam has an intensity profile as shown 
in Fig.^c), with a maximum intensity and thus barrier 
energy aligned with the center of the TOP trap. The 
barrier's potential energy decreases to zero over ~35 /im 
along three arms separated by azimuthal angles of ap- 
proximately 120°. With 170 [iW in the beam, which 
corresponds to a maximum barrier height of fcsx26 nK 
for our beam, three condensates are created and do not 
merge together during their growth; a set of three BECs 
created under these conditions is shown in Fig. QJd). 
With instead 45 fiW in the beam, corresponding to a 
maximum barrier energy of /csx7 nK, three independent 
condensates again initially form, but as the condensates 
grow in atom number, their repulsive interatomic inter- 
actions eventually provide enough energy for the three 
condensates to flow over the barrier arms. The three 
condensates thus naturally merge together into one BEC 
during evaporative cooling, as shown in Fig.^e). 

In our first study, three spatially isolated condensates 
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were created in the presence of a strong barrier of maxi- 
mum potential energy kgx26 nK, and were then merged 
together by ramping down the strength of the barrier to 
zero linearly over a time r. Any vortex cores formed by 
this process in the resulting BEC have a size below our 
optical resolution limit, and are too small to be directly 
observed in the trapped BEC. We thus suddenly removed 
the trapping potential and observed the atom cloud us- 
ing absorption imaging along the vertical axis after 56 ms 
of ballistic expansion. This entire process was repeated 
between 5 and 11 times for each of 6 different values of 
the barrier ramp-down time r between 50 ms and 5 s. 

In a significant fraction of our experimental runs, we 
observed one or more vortex cores in our condensates, a 
clear indication that condensate merging can indeed in- 
duce vortex formation. Moreover, the observed spatial 
density distributions vary from shot to shot, as would be 
expected with indeterminate phase differences between 
the initial condensates. However, many images are ab- 
sent of any vortices. Example images of expanded BECs 
in Fig. Ela)-(d) show the presence of vortex cores after 
various barrier ramp-down times. An analysis of vortex 
observation statistics is given in Fig. |2{e) for the differ- 
ent values of r examined. We define a vortex observation 
fraction F v as the fraction of images,for each value of t, 
that show at least one vortex core |2fj |. The error bars 
reflect the ambiguity in our ability to determine whether 
or not an image shows at least one vortex. For exam- 
ple, core-like features at the edge of the BEC image, or 
core-like features obscured by imaging noise, may lead 
to uncertainty in our counting statistics and determina- 
tion of F v . As the plot shows, F v reaches a maximum 
near 0.6 for the smaller r values, and drops to ^0.25 for 
long ramp-down times. We expect that with larger sam- 
ple sizes, F v should approximate P v for each r. Thus 
our results are consistent with our conceptual analysis, 
where P v > 0.25 for fast merging times when interference 
fringes may occur, and P v — 0.25 for long merging times. 

For r < 1 s, multiple vortices were often observed in 
our images, as images of Fig. [5] show, indicating that 
phase gradients are likely to play an important role in 
vortex formation if condensates are quickly merged to- 
gether. Furthermore, observations of multiple vortex 
cores may indicate the presence of vortices and antivor- 
tices. Although we are unable to determine the direction 
of fluid circulation around our observed vortex cores, we 
performed a test in which the barrier was ramped off in 
200 ms, thus forming multiple vortex cores with a high 
probability. We then inserted additional time to hold 
the final BEC in the unperturbed harmonic trap before 
our expansion imaging step. After such a sequence, the 
probability of observing multiple vortices dropped dra- 
matically: for no extra hold time, we observed an aver- 
age of 2.1 vortex cores per image, whereas this number 
dropped to 0.7 for an extra 100 ms hold time, suggestive 
of either vortex-antivortex combination or other dynam- 
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FIG. 2: (a)-(d) Images showing vortices in condensates cre- 
ated in the presence of a strong (fcsx26 nK) barrier. Prior 
to release from the trap, the barrier was ramped off over the 
time r indicated. Each absorption image shows a region of 
170 /im per side, (e) Vortex observation fraction F v vs. bar- 
rier ramp-down time r. The data for r values of 50 ms, 100 
ms, 200 ms, 500 ms, 1 s, 3 s, and 5 s, consisted of 5, 11, 10, 
10, 10, 5, and 5 images, respectively. For clarity, statistical 
uncertainties due to finite sample sizes are not reflected in this 
plot, but they generally exceed our counting uncertainties. 

ical processes by which vortices leave the BEC. However, 
images showing single vortices were observed even after 
5 s of extra hold time in our trap following the 200 ms 
barrier ramp, indicating a relatively long vortex lifetime 
in our trap. 

In a second investigation, we differed from the above 
experiment by using a weaker barrier with a maximum 
energy of /csx7 nK such that the three condensates nat- 
urally merged together into one BEC during the evapo- 
rative cooling process. We emphasize that this merging 
process is due solely to the increasing chemical poten- 
tials exceeding the potential energy of the barrier arms 
between the condensates; the barrier strength remained 
constant throughout the growth and merging of the con- 
densates when vortices may form. After evaporative cool- 
ing, we ramped off the optical barrier over 100 ms and 
released the atoms from the trap to observe the BEC af- 
ter ballistic expansion. Under these conditions, our vor- 
tex observation fraction was F lJ =0.56±0.06 in a set of 16 
images, with example images shown in Fig.^a) and (b). 

By adding an extra 500 ms of hold time after the fi- 
nal stage of BEC formation but before the start of the 
100 ms barrier ramp-down and ballistic expansion, the 
vortex observation fraction decreased to F„=0.28±0.14. 
Again, this drop in probability may be due to vortex- 
antivortex combination during extra hold time in the 
weakly perturbed harmonic trap. We thus conclude that 
with a maximum barrier energy of ks x7 nK, vortices are 
formed during the BEC creation process rather than dur- 
ing the ramp down of the weak barrier, consistent with 
our phase-contrast images of trapped BECs that show a 
ring- like rather than segmented final density distribution. 

Barrier strengths between the two limits so far de- 
scribed also lead to vortex formation, either during BEC 



4 



(a) 


<b) 


(c) 


(d) 

i 


□ 


| 


(g) 

i 
v. 


(h) 

• 


(i> 

* 


G) 

i 



FIG. 3: (a),(b) Images showing the presence of vortices natu- 
rally occurring in condensates created in a potential well with 
a fesx7 nK segmenting barrier, (c)-(h) Images obtained us- 
ing intermediate barrier heights. (i),(j) Images taken in the 
absence of a central optical barrier. Note the presence of a 
spontaneously formed vortex in (j). All absorption images 
represent a region of 170 fim per side. 



growth or during barrier ramp down. With a barrier 
strength in this range, up to at least four clearly defined 
vortex cores have been observed in single images, as the 
examples of Fig. EJc)-(f) show. Density defects other 
than clear vortex cores have also appeared in our im- 
ages, as in the upper left of Fig.^g), where a "gash"-like 
feature may be a possible indicator of vortex-antivortex 
combination; similar features have been seen in related 
numerical simulations |18| . Often, however, no vortices 
are observed; an example with no vortex cores appearing 
is Fig. 0th) . For comparison, an expansion image taken 
after creating a condensate in a trap without a barrier is 
shown in Fig.[3ji). 

Perhaps surprisingly, single vortex cores have also ap- 
peared in ~10% of our expansion images taken in the 
absence of any central barrier. In other words, for our 
basic single BEC creation procedure as outlined above, 
and without a segmenting barrier ever turned on, vor- 
tices occasionally form spontaneously and are observable 
in expansion images. An example image is shown in Fig. 
OUj). These observations may be related to predictions 
of spontaneous vortex formation due to cooling a gas 
through the BEC transition [2l|. We are currently in- 
vestigating these intriguing observations further. 

We finally note that to generate vortices by the mecha- 
nism described in this paper, it is important for two rea- 
sons that condensates merge and interfere while trapped, 
as opposed to during expansion. First, in a trapped BEC, 
the nonlinear dynamics due to interatomic interactions 
would play a key role in the structural decay of interfer- 
ence fringes, which may be responsible for generation of 
multiple vortices and antivortices seen with fast merging 
times. In an expanding gas, the interactions become neg- 
ligible as the gas density decreases. Second, by keeping 
condensates trapped during their mixing, we are able to 
study slow merging, where we believe that relative over- 
all phases are primarily responsible for vortex generation. 
We conjecture that with slow merging, it may be possi- 
ble to directly imprint relative phases onto three or more 
trapped, separated, and phase-correlated condensates to 



controllably engineer vortex states. 

In summary, we have generated vortices by merging 
together isolated and initially uncorrelated condensates 
into one final BEC. We have shown that our vortex ob- 
servation statistics are consistent with a simple concep- 
tual theory regarding indeterminate phase differences be- 
tween the initial condensates; however, quantitative the- 
oretical examination is needed for further analysis of our 
results. We have also demonstrated that condensates cre- 
ated in the presence of weak trapping potential defects 
or perturbations, such as our weak optical barrier, may 
naturally acquire vorticity and nonzero orbital angular 
momentum. This result challenges the common notion 
that a BEC necessarily forms with no angular momen- 
tum in the lowest energy state of a trapping potential; 
rather, the shape of a static confining potential may be 
sufficient to induce vortex formation during BEC growth, 
a concept that may be of relevance to other superfluid 
systems as well. 

We thank Ewan Wright and Poul Jessen for helpful dis- 
cussions, and Tom Milster for use of his Maskless Lithog- 
raphy Tool to create our optical barrier mask and phase 
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